Spin Glass Phase Transition on Scale-Free Networks 

D.-H. Kim,' G. J. Rodgers,^ B. Kahng,' and D. Kim' 

'School of Physics and Center for Theoretical Physics, 
Seoul National University, Seoul 151-747, Korea 
^Department of Mathematical Sciences, Brunei University, 
Uxbridge, Middlesex, UBS 3PH, United Kingdom 
(Dated: February 2, 2008) 

Abstract 

We study the Ising spin glass model on scale-free networks generated by the static model using the 
replica method. Based on the replica-symmetric solution, we derive the phase diagram consisting of the 
paramagnetic (P), ferromagnetic (F), and spin glass (SG) phases as well as the Almeida-Thouless line as 
functions of the degree exponent A, the mean degree K, and the fraction of ferromagnetic interactions r. To 
reflect the inhomogeneity of vertices, we modify the magnetization m and the spin glass order parameter 
q with vertex-weights. The transition temperature Tc (Tg) between the P-F (P-SG) phases and the critical 
behaviors of the order parameters are found analytically. When 2 < A < 3, Tc and Tg are infinite, and the 
system is in the F phase or the mixed phase for r > 1/2, while it is in the SG phase at r = 1/2. m and q 
decay as power-laws with increasing temperature with different A-dependent exponents. When A > 3, the 
Tc and Tg are finite and related to the percolation threshold. The critical exponents associated with m and q 
depend on A for 3 < A < 5 (3 < A < 4) at the P-F (P-SG) boundary. 

PACS numbers: 89.75.Hc, 75.10.Nr 
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I. INTRODUCTION 



Recently, considerable effort has been devoted to understanding complex systems by means 

of networks mm 

. An emerging phenomenon in real-world complex networks is a scale- 
free (SF) behavior in the degree distribution, Pdik) ~ k^^, where the degree k is the number of 
edges connected to a given vertex and A is the degree exponent 1^. Due to the heterogeneity of 
degree, many physical problems on SF networks exhibit distinct features from those in Euclidean 
space. For example, the critical behavior of the ferromagnetic Isingmodel on SF networks exhibits 
an anomalous behavior depending on the degree exponent A I^Q, 0, 0, While the critical 
behaviors are of the mean field type for A > 5, they exhibit an anomalous scaling for 3 < A < 
5. Moreover, the magnetization, m, decreases with increasing temperature as m ~ t-V(3-a) 
for 2 < A < 3, and so on fl, 0]. The Ising spin system on the complex networks, besides 
being of theoretical interest, can be used to describe various real world phenomena. For example, 
the two Ising spin states may represent two different opinions in a society. Depending on the 
interaction strength between neighbors, the overall system can be in a single or mixed opinion 
states, corresponding to the ferromagnetic or paramagnetic phase, respectively. 

In complex systems, such a description with only ferromagnetic interactions may not be suf- 
ficient in certain circumstances. In social systems, for example, the relationship between two 
individuals can be friendly or unfriendly. In biological systems, two genes can respond to an ex- 
ternal perturbation coherently or incoherently in microarray assay. For such cases, the spin glass 
model is then more relevant to account for such competing interactions. Recently, the spin glass 
problem has been studied on the small world network proposed by Watts and Strogatz llli through 
both the replica method and the cavity method Since SF networks are ubiquitous in nature, 
here we study the spin glass model on SF networks. 

The spin glass p roblem in the Euclidean space has been studied for a long time by various 
methods I3>13>[3>[3]. Most of the studies for spin glasses have concentrated on regular lattices 
or the infinite-range interaction model on fully-connected graphs, for example, the Sherrington- 
Kirkpatrick (SK) model Iryj . To achieve our goal here, we follow the study of the dilute I sing 



smng 
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ass model with infinite-range interactions, first performed by Viana and Bray (VB) 



lie i sing 

fly 



, because the model is equivalent to the Ising spin glass problem on the random 

nri 

graph proposed by Erdos and Renyi (ER) 112411251] . The ER random graph may be constructed as 
follows. The number of vertices N is fixed and assumed to be sufficiently large. Each vertex i 
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(i = 1, 2, . . . , A^) is assigned a weight pi, which is given as pi = 1/N, independent of the index i 
for the ER model. Two vertices i and j are selected with probabilities pi and pj, respectively, and 
if i 7^ j, they are connected with an edge unless the pair is already connected, which we call the 
fermionic constraint. This process is repeated NK/2 times. In such networks, the probability that 
a given pair of vertices {i, j) (i ^ j) is not connected by an edge, denoted by 1 — fij, is given by 
(1 — 2piPj)^^/'^ ~ ex-p{—NKpiPj), while the connection probability is 

f,, = l-exp{-NKp,p,). (1) 

Since PiPj = for the ER graph, the fraction of bonds present becomes fij ^ K/N and the 

average number of connected edges is NK/2. So K is the mean degree, and corresponds to p of 
Ref. 1 18]. 

The SF network can be constructed through a generalization of the above to the case where the 
vertex- weights are given by 

Pi = -^T^i (2) 

where is a control parameter in the range [0, 1), and CAf(/^) = Yl!j=ij~^ ~ N"^^^ /{\ — /i). 
Then the resulting network is a SF network with a power-law degree distribution, Pd{k) ~ k~^, 
with A = 1 + The model is called the static model, where the name 'static' originates 

from the fact that the number of vertices is fixed from the beginning This model has the 
advantage that many of its theoretical quantities can be calculated analytically H. Note that 
since NKpiPj ~ N^''-^/{ijY for finite K, when < ;U < 1/2 (A > 3), 

^ NKp,p^, (3) 

however, when 1/2 < < 1 (2 < A < 3), /jj does not necessarily take the form of Eq.©, but it 
is given as 

{1 when ij < A^2-i/m^ 

(4) 
NKpiPj when ij > iV^"^/^. 

This is due to the fermionic constraint that at most one edge can be attached to a given pair of 
vertices. The mean degree of a vertex i is NKpi and the mean degree of the network is K 11271] . 

In this work, we study the Ising spin glass model defined on the static model. In Sec. II, we in- 
troduce the Hamiltonian of the spin glass system on the static model and derive the free energy by 
using the replica method. We also introduce physical quantities such as the magnetization, and the 



spin glass order parameters in a modified form. In Sec.III, we present the replica-symmetric solu- 
tions by using the SK-type approximation, from which the phase diagram including the Almeida- 
Thouless line and the critical behavior of the spin glass order parameters are derived. In Sec. IV, 
we use the perturbative approach to derive the phase diagram and the critical behaviors of the or- 
der parameters, and compare them with those obtained from the SK method. The final section is 
devoted to the conclusions and discussion. 



n. THE SPIN GLASS MODEL 



We consider the Ising-type Hamiltonian, 

= - ^ JijSiSj {si = ±1), (5) 

defined on a graph G realized by the static model. Jij is nonzero only when the vertices i and j 
are connected in G. The network ensemble average for a given physical quantity A is taken as 

{A)k = J2Pk{G)A{G), (6) 

G 

where Pk{G) is the probability of G in the ensemble and (• ■ ■)k the average over different graph 
configurations. For the static model we consider here, it is given that 

pk{g) = n k n (1 - k) (7) 

with fij = 1 — exp{—NKpiPj), pi being given in Eq.©. 

In the spin glass problem, the coupling strengths {Jij} are also quenched random variables. We 
assume in this paper that each J^j is given as + J or — J with probability r and 1 — r, respectively, 
so that the coupling strength distribution is given as 



(8) 



Pr{{J^J})= n [rS{J^J-J) + {l-r)6{JiJ + J) 

The case of r = 1 (r = 1/2) is pure ferromagnetic (fully frustrated) one, and we consider r in the 
range ofl/2<r< 1 throughout this work. The average of a quantity A with respect to Pr{{Jij}) 
is denoted as {A)r. Thus the free energy is evaluated as —^F = {{InZ)^)^ with Z being the 
partition function for a given distribution of {Jij} on a particular graph G. 
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In this paper, the replica method is used to evaluate the free energy, i.e., — 
lim„^o[((^")r)i^ — To proceed, we evaluate the n-ih. power of the partition function. 



{i..j)<^G «=1 

n 

= Tr{^c,}exp 



i<j a=l 

n 

'^^•)-l)J}], (9) 



where the trace Trjgc,} is taken over all replicated spins s" = ±1, a = l,...,nis the replica 
index, and /? = 1/T. We mention that the disorder averages over Pk{G) and Pj.{{Jij]) can be 
done simultaneously since both types of disorders are independently assigned to each edge of the 
fully-connected graph of order A^. The part inside the exponential in Eq.® can be written in the 
form, 

n n 

5^ In {l + (( exp(/3J,,. 4s]) - l) J } = E NKp,p,(^exp{(3J,, " + ^ 

i<j a=l i<j a=l 

(10) 

where R stands for the remainder which are of higher order in K. It is shown in APPENDIX 
A that for finite K, Eq.^ is 0{N) while R is at most C(l) for A > 3 and 0{N^^^ InN) for 
2 < A < 3, so that it can be neglected in the free energy calculation. 

Once R in Eq.lfTOt can be neglected, we can proceed as in VB lisl . By using the relation, 

n 

( exp(/5J,, J2 ))^ = ( n [cosh(/5J,,)(l + s^sj tanh(/5J,,))] (11) 

a=l a 

in Eq. lfTOb and applying the Hubbard-Stratonovich identity, Eq.® is reduced to the form 

{{Z^)r)K = j rfqexp{-iVn/?/(q)}. (12) 
The intensive free energy /{q}(= F/Nn) in the thermodynamic limit (A^ oo) then becomes 
nmq} = ^E^' + ^E^'/3 + ^ E ?'/37 + ----^El^Tr|,.|expX„ 

a a<f3 Q!</3<7 i 

(13) 

where 

X, = NKT,p,Y,1aS? + NKT2P^J2qapsts^ + NKTsP, Yl ^c^m^IsUJ + ' ' ' ^ (14) 

a a<P a</3<7 



and 



Ti{T) = (cosh"/5Jijtanh'/5J,j)^ ^ [r + (-l)'(l - r)] tanh' /5J (/ = 1,2,...)- (15) 

TrjsQ} is the trace over the replicated spins at vertex i and the N ^ oo limit is to be implicitly 
understood to the expression The elements of a set {q}, g^,, g^^, qa/s-y, etc., defined as 

(la = ^Pi{st)i, Qap = ^Pi{SiS^)i, Qaf^^ = ^ Piisf s]) i, CtC. (16) 

i i i 

are the order parameters of the spin glass system, called the magnetization, the spin glass order 
parameter, and so on. The average is evaluated through {A)i = Tr{5a|AexpXj/Tr{sa| expXj. 
Note that unlike the case of the ER random graph, the order parameters are summed with weight 
{pi} in Eq. lfT^ due to the inhomogeneity of the SF networks. For the ER case however, pi = 1/N 
and it becomes that g„ = Y.i{sf)i/N, q^p = Y.i{s?s^)i/N, q^p^ = Y.i{s?s^ s])i/N, and so on 
llisll . To distinguish, we use bar notation for the unweighted cases. 

Here we consider the replica symmetry (RS) in which spins with different replica index are 
indistinguishable, and we invoke two methods to determine the phase boundaries of the ferromag- 
netic (F), paramagnetic (P) and spin glass (SG) phases and the temperature dependences of the 
order parameters. The first is the approach similar in spirit to SK in which higher-order terms than 
qai3 in Eqs. (fT3b and (fT4b are neglected. In this method, the remaining two order parameters as well 
as the Almeida-Thouless line can be obtained for all temperatures. The second is the perturbative 
approach used in VB. In this case, we expand the term of InTrexpXj in Eq. (fT3b up to appro- 
priate orders, and the order parameters qa, q^p, qap-y and qap-yS are explicitly calculated. Through 
the perturbative approach, we can find that the contributions by higher order terms such as qap-y 
are negligible compared with those by qa and qap near the phase transition points. Thus, the two 
methods produce identical results for the phase boundaries and the same critical behaviors near 
the transition points for the two order parameters, qa and qap. 

III. THE SHERRINGTON-KIRKPATRICK APPROACH 
A. The replica symmetric free energy 

We first study the RS solution IitI and obtain the phase boundaries of P, F and SG. For sim- 
plicity, the RS magnetization and the RS spin glass order parameter are denoted as m(= g^) and 
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q(= Qa/s), respectively, and the free energy expression Eq. (fT3t is truncated at the order of q. Then 
the RS free energy is rewritten as 

, KTi 2 KT2n(n-l) ^ 1 ^ 

i 

with 

Z, = Tr|,.} exp \^NKT,p,m J^^? + NKT^p^q^^^^^^-^]. (18) 

a 

By using the Hubbard-Stratonovich identity, Zi can be rewritten as 

Z, = exp{--NKT2P^q} / Vz[2 cosh r]i{z)r, (19) 

where / Vz ■ ■ ■ = dz e~^^/^ • ■ ■ and rii{z) = NKTiPim + z\/NTCT^^p~q. Then in the 

limit of n — > 0, the RS free energy becomes 

A'" 

(3f{m,q) = ]^KT,m^ + ]^KT2q--^KT2q^ - j ^ ln[2 coshr/,(^)]. (20) 

i=i 

By applying {df /dm) = and (df/dq) = to the free energy, Eq.(l20b. we obtain the coupled 
self-consistent equations for m and q to be 

AT 



m 



I Vz^Pitsxih{NKTiPim + z^/NKT^), (21) 

i=l 



and 

N 



q= Vz^pi tauh\NKTiPim + zs/NKT^). (22) 
^ i=i 

In Eqs. (|2Tl) and (l22b . we can see that q cannot be zero unless both m and q are zero, while m can 
be zero even when g 7^ which defines the SG phase. 



B. The phase boundaries 

The P-F (P-SG) phase boundary is given as the temperature, the Curie temperature (the spin 
glass phase transition temperature Tg), where m (g) starts to be nonzero. We first consider the case 
of A > 3. When m and q are small, the free energy, Eq. (l20b . is written as 

i i 

+ higher order terms. (23) 
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FIG. 1: The phase diagram in the {K/Kp,T/J) plane for A = 4.5(> 3) with r = 1/2 (a) and A = 4.5(> 3) 
with r = 2/3 (b), and the same in the {K,T/J) plane for A < 3 with r = 1/2 (c) and A = 2.5(< 3) with 
r = 2/3 (d). Note that A'p = for 2 < A < 3 in the thermodynamic Umit. 

It is known that as K increases, the static model undergoes the percolation transition at 

1 (A-l)(A-3) 



(24) 



Since ^ . p: 



2 



NE^Pl (A -2)2 • 
{{k'^)K — {k)K)/{k)\ with {k)K = K and (fc^)x denoting the first and the 



second moments of the degree for a give n mean degree K, respectively, Eq. (|24|l is equivalent to 
the condition {k'^)K = '^{k)K 



;e for a give r 



Ti(r,) 



Thus one obtains that 

Kp/K for P-F, and 
Kp/K for P-SG, 



(25) 
(26) 



where Ti(r) = (2r - 1) tanh( J/T) and T2(T) = tanh^{J/T). Note that when K/Kp < 1, 
there is no solution of Eqs. (l25b and (l26b . implying that the system is always in the P state. This is 
because the network has an infinite component only for K > Kp. When r = 1/2, Ti = and the 
phase diagram is rather simple. The P-F transition does not occur, and the system is either in P or 
SG phase whose boundary is given by Eq. (l26t . FIGda) is the phase diagram in the {K/Kp, T / J) 
plane for the fully frustrated case (r = 1/2) for A > 3. When 1/2 < r < 1, both the F phase 
(m 7^ 0, g 7^ 0) and the SG phase (m = 0, g 7^ 0) appear. FIGHIb) is the phase diagram for a 
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partially frustrated case with r = 2/?) and A = 4.5, which is a prototypical case ofl/2<r<l and 
A > 3. For K/ Kp < 1, only the P phase appears, but for K/ Kp > 1, several phases exist. There 
exists a multicritical point {K* /Kp, T* / J), where the P-SG-F phases merge, which is determined 
to be 

(27) 



- 



\Kp' J J V(2r-l)2'tanh-^(2r-l^ 
by setting Ti(T*) = T2(T*) = Kp/K*. For Kp < K < K*,theP phase goes into the SG phase, 

while it goes into the F phase for K > K* as temperature is lowered. As r ^ 1, the multicritical 
point converges to (1, 0), indicating that only the P-F phase transition occurs. As r — > 1/2, it shifts 
to (oo, oo), indicating that only the P-SG phase transition occurs as shown in FIGlUa). 

Besides the P, F, and SG phases, the mixed (M) phase is present, which is defined as the re- 
entrant SG phase with nonzero macroscopic ferromagnetic order, located below the F phase as 

HQ 

temperature is lowered II15LI16I1 . The SG-M phase boundary is determined as the vertical straight 
line from the multicritical point to T/ J = [.301 . The F-M phase boundary is determined by the 
so-called Almeida- Thouless (AT) line IbiI . 

r ^ , 

{KT2y' = rDzJ^ scch\NKT^pim + z^/NKT^), (28) 
i=i 

which is obtained easily by multiplying vertex- weights to the AT line formula of the SK model, m 
and q above are the solutions of Eqs. lOn) and (E^ . We determine T satisfying Eq.(EHt numerically. 
The F-M boundary in FIGHIb) exhibits a fat-tail behavior, implying that the M phase persists for 
large K. This AT line is the phase boundary between the replica symmetric phase and the replica- 
symmetry-broken one. Thus, Eq. (l28b indicates the region where the replica-symmetric solution 
derived in the following sections is valid. We also check the P-SG boundary from Eq.(EHl). which 
is the same as Eq.lE^. 

Next we consider the case 2 < A < 3. In this range, Kp ~ 7V-(3-^)/{-^-i) ^ g as ^ oo 
and consequently Tc and Tg ^ oo. Thus the whole {K, T / J) plane is covered with the ordered 
states. FIGHJc) is the phase diagram for the fully frustrated case (r = 1/2) for A < 3. The P phase 
appears only for = 0, and the SG phase is located in the region > 0. FIGHId) deals with the 
case of 1/2 < r < 1 and A < 3. The P phase appears only at K = 0, but for > the F and 
M phases appear and the F-M boundary is given by the AT line (Eq.(ESt). As r ^ 1, the M phase 
disappears and only the F phase appears in the region of > 0. 

We also consider the phase diagram in the (r, T / J) plane for given A and K in FIGEI The 
phase diagram is schematically similar to the one for the SK model. In the original paper of 
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FIG. 2: The phase diagram in the (r, T/ J) plane for = 5 with A = 4.5(> 3.0) (a) and A = 2.5(< 3.0) 
(b). 

the SK model B, a new coupling constant Jq of the F interactions was introduced and the ra- 
tio JqI J plays a similar role of the parameter r here. Accordingly, the phase diagram in the 
(r, T I J) plane here corresponds to the one in the (Jq/ J, T / J) plane in the work of the SK model. 
FIG|3a) shows the phase diagram for A > 3. The formulae of the phase boundaries of P-SG 
and P-F are easily derived from Eqs.(E5l and (E^ . The P-SG phase boundary is constant as 
1 / tanh"^ a/ Kpl K, independent of the parameter r and the P-F phase boundary is determined as 
T / J =1/ taii\i~^[Kp/ K{2r — 1)). The multicritical point is determined as 

2 'tanh-,/^> ^^^^ 
The SG-M phase boundary is given by the vertical line as before. The F-M boundary is obtained 



from Eg. (1281) . finding numerically that the region of the M phase shrinks as A increases, and 
eventually it remains on the line spanning from the multicritical point to T = for a given K, 
while it exhibits a fat-tail behavior in the direction of the parameter K. 

We plot the phase diagram in the (r, T / J) plane for A < 3 with a given K{> 1) in FIGlZtb). 
Note that as A ^ 3 for a given K, r* approaches 1/2, while T* / J diverges to infinity. Thus, for 
2 < A < 3, the SG phase can exist only when r = 1/2. For 1/2 < r < 1, the F and M phases exist 
and the F-M boundary is given by the AT line (Eq.dSSt). 



C. The SG order parameter 

In the SG phase (m = 0, g 7^ 0), the SG order parameter q is determined by 



N 

q= / p2^pitanh^(;zv^iV^T^). 



(30) 
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Note that Eq.dSOb is independent of r but valid for 1/2 < r < r*, r* being the value of r at the 
multicritical point. 

In this section, we determine the critical behavior of q near the SG transition. The right hand 
side of Eg. (BUI) involves a sum of the type 

1 ^ 

1=1 

with y = {1 — fi)KT2qz'^ and F{x) = xtanh^ y/x. When y is small in S{y), a singular term 
y^^^ competes with other regular terms. General expressions for small y expansions are derived in 
APPENDIX B. When Eq. (IB7l) is used and the Gaussian integration over z is performed, Eg. (BUI) 
becomes 

q/{l - /i) = ^2^-ir(A - l)D{X)Q"-' - ^Q' + 2^Q" + 0{Q'') (32) 

where 

f dx tanh^ x for 2 < A < 3, 

D{X) = { ■''^ (33) 

I - dx x3-2^(a;2 - tanh^ x) for 3 < A < 4, 

and Q' = (1 - fi)KT2q = (A - 2)irT2g/(A - 1). Equating the right hand side of Eq.lEl) with 
g7(l - ^ifKT2, one sees that Q'^-^ ~ KT2 ~ for 2 < A < 3, [l/Kp - I/KT2) ~ 
for 3 < A < 4, and {l/Kp- l/isTTa) ~ Q' for A > 4. Here fsTp is given by Eq.dSll) and 
the A-dependent positive coefficients are neglected. Therefore, as T ^ 00 (2 < A < 3) or 

Eg = {Tg - T)/Tg ^ (A > 3), ^ bchavcs as 

^ y-2(A-2)/(3~A) fo^ 2<A<3, 

ey^^-') for 3<A<4, (34) 

eg for A > 4. 

When A = 3, use of Eq.dHHl) yields 

g ~ exp{-2TyKJ^) as T ^ oo, (35) 

while, when A = 4, 

g~eg/lneg^ as — >• 0. (36) 

For general temperatures, q can be obtained numerically from Eq. (l30b . The behavior of q for 
various A are shown in FIGjSl 
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FIG. 3: The behavior of q in Eq.(|30|l for N = 5000 and = 2 for 1/2 < r < r* 
IV. THE PERTURBATIVE APPROACH 



In this section, we use the perturbative approach to evaluate the free energy and to obtain 
the order parameter behaviors near the transitions. For simplicity, we use the notations defined 
through Qc, = KT^q^, Q^fi = KT2qap, Qap-y = KTsq^p^, Q^fj^s = KT^q^p^s, and so on. 
Let R represent a subset of the replica indices {1, 2, . . . , n}. Then it is convenient to denote the 
set {Qa, Qaf3, . . .} as {Qr}. We also write ctr = HaeR ^ With these notations, Eq. (fT4l) 
becomes Xi = NpiQ^a^ where the sum is over all subsets of {1,2,..., n} except the null 
set, and 

= J] e^^''««'^« = H cosh Np.Qr J](1 + T^a^) (37) 

R, R R, 

with tr = tanh NpiQ^. Our perturbative approach is to expand Ylni^ + trctr) and keep only the 

1/2 

term s up to given order. In the ER limit A ^ oo, we anticipate that ~ Cc , Tq,^ ~ ec, etc. from 
VB lIlSll . where ec = (Tc — T)/Tf. is the reduced temperature. 

Using the properties that Tr ctr = 0, TrcTRCXR/ = for R 7^ R' and so on, the first few terms 
relevant to our discussion below are 



a a<l3 q:</3<7 q</3<7<<5 

i R i a</3 Q</3<7 

a</3<7 a</3<7 
a<P<-y<5 a</3<7<(5 



2 

0/375 
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a<P<-y<5 



The result of APPENDIX B with F{x) = In cosh x gives 
i 5^1ncosh(iVp,QK) = A(A)Q^-^ + ^Ql 

i 

where 



A{X) ^ { 



(A- 


1)A 


-2 


(A- 


2)^ 


-1 


(A- 


1)^ 


-2 


(A- 


2)^ 


-1 



L (A-l)^ 



Jg°° (ia; s ^ In cosh x 

dx x"^ ^ In cosh x 
Jq°° (ia; ^ In cosh x 



1™2 
2^ 



ir2 -L J-r" 
2 12 



for 2 < A < 3, 
for 3 < A < 5, 
for 5 < A < 7, 



and 



az = (A-2)7[(A-l)'-i(A-l-0]. 
The last sums in Eq. (l38l) can be represented as integrals as 

— trTr' ■ ■ ■ = (A — 1) / dz tanh zQ'-^ tanh 

withQ'j.^ (A-2)Qr/(A-1). 



(38) 



(39) 



(40) 



(41) 



(42) 



A. The replica symmetric free energy 

We derive the RS solution of the order parameters up to the fourth order with the notations of 
Qa = M, Qaf3 = Q, Qafs-y = Qs and Qap-yS = Qa, respectively. Then the terms in Eq.lH^ take 
the form of 

/oo 
dz z-^ tanh'^i zM' tanh''^ zQ' tanh"^ zQ'^ tanh"^ zQ'^, (43) 

where M' = (1 — /u)M = (A — 2)M/ (A — 1), ni, . . . , n4 are integers, and other primed quantities 
are similarly defined. 

The RS free energy f{M,Q,Q3,Q4j in the limit of n ^ is then written as 



6 



12 



24 



36 



48 



-A{X) 



4 



(44) 



+ 2^2,1,0,0 — -^[^a.oa.o + 3^2,2,0,0 + ^0,3,0,0] + ^[^4,0,0,1 + 3i3o,2,o,i + ^0,4,0,0] 
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where bi = {KTi)~^ — 02 for / = 1, 2, 3 and 4. Note that 1/02 is nothing but Kp for A > 3 given 
in Eq. (l24b . while it is negative for 2 < A < 3. 

The RS solutions of M, Q, and Q4 are obtained by solving the self-consistent equations, 

df/dM = df/dQ = df/dQs = df/dQ, = 0. (45) 

When M, Q, and Q4 are small, ... are small. Their leading order behaviors are calculated 
in APPENDIX C. 

The phase boundary of the P-F transition is determined as the same obtained in the SK ap- 
proach. When 2 < A < 3, since A{\) is nonzero and positive, the transition temperature 
becomes infinity so that the system is always in the F phase when r > 1/2. For r = 1/2, however, 
bi = 00, and has to be zero. Then the system is in the SF phase. 

B. The P-F transition and the order parameters 

We first consider the P-F transition. In the F phase, all the four order parameters remain 
nonzero. The behaviors of each order parameter within leading order are discussed below and 
listed in TABLE I. 



(i) When 2 < A < 3, the leading order terms in free energy f3f are 

Pf ~ -A(A)M"-i + + (1a(A) - C2,o - ^C3,o + ^C4,o + ■ • • )q 



order parameters 


2 < A < 3 


3 < A < 4 


4 < A < 5 


5 < A < 6 


A > 6 


M 


^ 2^-1/(3- A) 






1/2 
~ ec' 


1/2 
~ ec 


m 


^ 2.-(A-2)/(3-A) 


l/{A-3) 


l/(A-3) 
ec 


1/2 
~ ec 


1/2 
~ ec 


Q 


^ T~(4-A)/(3-A) 


(A-2)/(A-3) 

~ ec 


2/(A-3) 
ec 


~ei 


-e^ 


q 


^ 2^-(A-2)/(3-A) 


(A-2)/(A-3) 

~ ec 


2/(A-3) 
~ ec 


~ei 


~e^ 


Q3 


^ 2.-(7-2A)/(3-A) 


(A-2)/(A-3) 

~ ec 


(A-2)/(A-3) 

~ ec 


3/2 

~ ec 


3/2 

~ ec 


Q3 


^ 2.-(A-2)/(3-A) 


(A-2)/(A-3) 

~ ec 


(A-2)/(A-3) 

~ ec 


3/2 
~ ec 


3/2 
~ ec' 


Qa 


^ 2^-(10-3A)/(3-A) 


(A-2)/(A-3) 


(A-2)/(A--3) 

~ ec 


~ ,{A-2)/2 


~ e2 


Q4 


^ 2^-(A-2)/(3-A) 


(A-2)/(A-3) 


(A-2)/(A-3) 

~ ec 


~ ,{A~2)/2 


~e2 



TABLE I: The A-dependent critical behaviors of the four order parameters and their scaled quantities 
(Eq.dT/l) under the P-F transition. Here ec = (Tc — T)/Tc is the reduced temperature. 
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+^C2,iM^-2g - ic3,iM^-2Q3 + ic4,iM^-2Q4 (46) 
from TABLE III with C„ p given in Eq. (IC6l ). 

By applying Eg .(1^31 to the free energy, we obtain the self-consistent equations for the four 
order parameters. Note that from the definition of bi = (KTi)^^ — 02, we find that bi ~ 
as T — > 00. All other coefficients such as A(\) and {Cn,p} are independent of T. From 
df/dM = 0, we obtain that -(A - 1)A{X)M^-^ + biM = 0, leading to that M ~ [(A - 
l)A(A) /6i]i/(3-A)] _ 2--i/(3-A)_ pj.qj^ gj- ^ 0, we obtain that (A(A) - 2C2,o + ■ • •) - 
+ C2,iM'^^^ — 62Q = 0. Since the second term is more dominant than the first, 
we obtain that Q ~ C2,iM^''^/b2 ~ T'{*-^)/{'i->-) _ Fortunately, the coefficient of Q^'^ is 
not needed to determine the leading order behavior of Q. Similarly, we obtain that Q3 ~ 
7--(7-2A)/(3-A) ^ 2.-(io-3A)/(3-A)_ Subsequently, we obtain m q q^, ^ 

y_(A-2)/(3-A)^^here 

m = M/KTr, q = Q/KT2, q^ = Q^/KT^, and q^ = Qi/KT^. (47) 

It is noteworthy that the behavior of m is different from that of the unweighted magnetiza- 
tion, fn ~ T^^/''^^^\ where m = (l/N) Xlii-^i) as previously studied in Ref.Jvllsl]. This is 
because m ~ M to the leading order. 

(ii) When A > 3, the transition temperature Tc is determined by 

6i(r,) = 0, i.e., a2KT,{T,) = l (48) 
which is the same as Eq.(E5t. When 3 < A < 4, the leading order terms in Pf are 

/3f ^ jM' - A{X)M'-^ - + ic2,iM^-2g 

+ ^iQl - \c,,iM'-'Q, - ^iQl + Wm'~'Q,. (49) 
D 3 8 4 

Note that A(\) < for 3 < A < 5. The most leading term is (6i/2)M^ and the transition 
temperature Tc is determined by bi = 0. Just below bi < and ~ 0(ec), where 
ec = (Tc - T)/Tc. From df/dM = 0, we obtain that -(A - 1)A{\)M^-^ + biM = 0, 
leading to that M ~ e^^^^^^ From df/dQ = 0, we obtain that C2,iM^-2 - b2Q = 0. 
Since 62 is constant near Tc, we obtain that Q ~ Af^"^ ~ g^(A-2)/(A-3) Similarly, it is 
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obtained that Q3 ~ M^'^ ~ q(a-2)/(a-3) ^ ^ g^(A-2)/(A-3)_ unlike the case 

of 2 < A < 3, m ~ M, g ~ Q, ?3 ~ Q3, and 54 ~ Q4. Such relations hold for all A > 3. 



(iii) When 4 < A < 5, the free energy is written as 

(3f ^ - A(X)M'-' - + ^M'Q 

+ ^jQl - lCs,iM'-'Qs - + lc,,,M'-'Q,. (50) 

Following the same step as used in 3 < A < 4, we obtain that M ~ e^'^'^^'^^ Q ~ ~ 
e'J^'-'^ and Qs ^ Q, ^ M'-' ^ ,(a-2)/(a-3)_ 

(iv) When 5 < A < 6, the free energy is written as 

+^Q^ - ^M^ga - + ]C,,,M'''Q,. (51) 
o 3 8 4 

Following the same steps as before, we obtain that M ~ e^^, Q ~ ec, Q3 ~ ec^^, and 

(v) When A > 6, the free energy is written as 

/3/ ^ ^M^ + ^M^ - + ^M'Q + - ^M'Qs - ^Ql + ^M'Q,. (52) 
A-j 2 12 4 2 6^3 ^-^8 4 ^ ^ 

Using the same step as before, it is obtained that M ~ e^^, Q ~ ec, ~ ec^^ and Q4 ^ e^. 

It is interesting to note that as A increases, the order parameters progressively acquire the classical 
mean field behavior Qn ~ e"'^^ starting from the lower order ones. 



C. The P-SG transition and the order parameters 

Here we consider the P-SG transition. In the SG phase, M and Q3 are always zero for all 
temperatures. Thus, the free energy becomes simpler compared with that in the F phase. Using 
the same method as used in the P-F transition, we obtain the P-SG transition temperature and the 
order parameters Q and Q4 in various region of A, which is listed in TABLE n. 

For more details, we first determine the P-SG phase boundary. When 2 < A < 3, since A(X), 
the coefficient of Q^^^ is nonzero for all T, the spin glass transition temperature Tg is infinity, and 
no P phase exists for all T. When A > 3, the transition point Tg is determined by the formula 

b2iTg) = 0, i.e., a2KT2{Tg) = 1, or KT^iTg) = (53) 
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which is the same as derived in the SK method. In the SG phase, the order parameter behaves as 
follows: 

(i) When 2 < A < 3, the leading order terms of(3f read off from TABLE III with M = Q3 = 
are 

Pf ^ (^A(A) - ^Cs,o + ^C4,o + ■ ■ ■)q'-' - 

+^ [a{\)Q\-^ - ^Ql + 3C2,iQ^-'Q,) . (54) 

By applying df/dQ = df/dQ^ = 0, we obtain that Q ~ r-2/(3-^) and Q4 ~ Q^-^T^ ~ 
7--(8-2A)/(3-A) Using the relation Q = KT2q and (^4 = KT^q^, we obtain that q ~ 
g4 ~ 2^"2(A-2)/(3-A) rpj^g result of q is the same as the one derived through the SK method, 
Eq. (l34b . Note that the coefficient of Q^^^ in the perturbative approach is in the form of 
infinite series while the same is obtained in a closed form in Eq. d^ . 

(ii) When 3 < A < 4, the free energy is 

/?/ ^ -|q2 ^ [^(^)/2 - C3,o/3 + C4,o/4]Q^-i - + Ic^^^q>^-^q^, (55) 

We note that the coefficient 62 ~ —^g with tg = {Tg — T)/Tg. Then we obtain Q ~ e^^^'^"'^^ 
Similarly, from df/dQ^^ = 0, we obtain Q4 ~ Q^^^ ~ g^->^^2)/(A-3) ^^^.j^ being constant. 

(iii) When A > 4, we have 

Pf - - jQ^ - - + la,Q'Q.. (56) 
By following the same step above, we obtain that Q ^ eg and Q4 ~ e^. 



order parameters 


2 < A < 3 


3 < A < 4 


A > 4 


Q 
Qi 

94 


^ 2^-2/(3- A) 
^ 2.-2(A-2)/(3-A) 
^ 2^-{8-2A)/(3-A) 
2^-2(2- A)/(3- A) 


l/{A-3) 

~ eg 

l/{A-3) 

~ eg 

(A-2)/(A-3) 

(A-2)/(A~3) 
~ ^9 


^9 
^9 
9 



TABLE II: The A-dependent behaviors of the two order parameters and their scaled quantities in Eq. (l47l 
under the P-SG transition, where eg = {Tg - T) /Tg. 
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V. CONCLUSIONS 



We have studied the spin glass phase transition on SF networks through the static model. The 
model contains generic vertex-weights in it, and edges between two vertices are connected with 
the probability given in Eqs. ([T]) and ©. The static model enables one to study the spin glass 
problem using the replica method by generalizing the dilute Ising spin glass model with infinite- 
range interactions. Here we obtained the replica-symmetric solutions through the two methods, the 
Sherrington-Kirkpatrick approach and the perturbative approach. We also found the phase diagram 
consisting of the paramagnetic (P), ferromagnetic (F), spin glass (SG), and mixed (M) phases in the 
space of temperature T, the mean degree K, the fraction of the ferromagnetic interactions r, and 
the degree exponent A. The AT line was also obtained numerically. The phase diagram is shown 
in the {K, T) and (r, T) planes, which are presented in Figs. 1 and 2, respectively. The critical 
temperatures and Tg for the P-F and P-SG phase transitions are simply related to the percolation 
threshold Kp in Eqs. (l25t and (l26b . We obtain the same results in the two approaches. Thus Tc and 
Tg are infinite when 2 < A < 3. The magnetization and the spin glass order parameter are modified 
to account for the inhomogeneity of vertex degrees as m = J2iPi{^?)i ^^'^ ? — '^iPi{^'i^i)i-' 
where Pi is the weight of vertex i. Such quantities depend on the degree exponent A. When 
2 < A < 3, due to the fact that Tc = oo and Tg = oo, m and q decay as power-laws for large 
T as shown in TABLES I and II, which is different from the patterns of rh and q, defined with 
Pi = 1/N. When A > 3, the order parameters exhibit continuous phase transitions across T^ 
and Tg, and the associated exponents depend on A, which are listed in TABLES I and II. As Q3, 
Q4, ... are of higher orders, the SK approach in Sec. Ill, and the perturbative one in Sec. IV give 
the identical results for m and q to the leading order. We find the critical exponents for the P/SG 
transition are non-classical in the range 3 < A < 4, corresponding to 3 < A < 5 for the P/F 
one J^i. We have not presented our results at integer values of A in Section IV for simplicity. At 
the borderline cases of A, the logarithmic corrections as given in Eqs. (IB8I) . (IC3I) and (IC7I) should 
be considered explicitly. We mention that the finite-size effect is an important issue especially for 
2 < A < 3 which we leave for a further study. 

It is noteworthy that the method we developed here can be applied to other problems in equi- 
librium statistical physics on SF networks. A novelty in this approach is that one needs not rely 
on the local treelike structure of SF networks used e.g. in [Q]. The result of the phase diagram 
and the behavior of the order parameters may be helpful in understanding emerging patterns in 
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various systems with competing interactions such as social or biological systems. For example, in 
the region 2 < A < 3 where most real-world SF networks belong, it is known that the structural 
characteristic of the network is so dominant that homogeneously interacting systems are in the 
ordered state for all temperatures. Our result shows that it is also the case even when there are 
competing interactions. Also for 2 < A < 3, the fact that a slight dominance of cooperative inter- 
actions (r > 1/2) drives the system into the ferromagnetically ordered or the mixed state suggests 
that most social and biological systems would be driven into the majority state (ferromagnetic or 
mixed state) at equilibrium. While the current study is meaningful as a first step of understanding 
thermodynamic property of the systems with competing interactions, further studies have to be 
followed towards real-world systems where the signs of interactions may be correlated with the 
degrees of vertices, or the interaction signs may change with time as in the prisoner's dilemma 
problem. 

While preparing this manuscript, we have learned of a recent preprint by Mooij and Kap- 
pen which addressed the same issue. They used the Bethe approximation to obtain a crite- 
rion for Tg and applied it to the A = 00 and A = 3 cases numerically. Our work gives analytic 
results for Tg as well as physical ones such as the phase diagram and the behaviors of the order 
parameters, which depend on the degree exponent. 
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APPENDIX A: EVALUATION OF THE REMAINDER 

In this APPENDIX A, we show that 



with R < 0{N^-^\nN) for 2 < A < 3, i? < 0{{\nNf) for A = 3 and R < C(l) for A > 3. 
Here Sij = {exp{(3Jij J2a=i ^i^'j) ~ is a quantity independent of the system size N. To do so. 




(Al) 
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we expand the logarithm on the left hand side of Eq. dAll) to write it as 

i<j i<j i<j n=2 i<j 

and show that the positive quantities defined by 

R' = I J2iNKp,pj - f,j)S,,\ (A3) 

i<j 

and 

Rn^\J2f-jSr,\ (A4) 

i<j 

(n > 2) are all bounded above by o{N) quantities. 

First let us consider R'. Since Sij are independent of N, we replace Sij by their maximum 
value Smax = maxi<j \Sij \ to get 

R' < 5^ax E GliNKp.p,) <^[J2 Gi{NKp,p,) - G,{NKpl)], (A5) 

i<j id 

where 

Gi{x) =x -l + e"". (A6) 

Here we have added i = j terms for i > 2 on the right hand side of Eq . (IAS I) for convenience. Since 
Gi{x) is monotone increasing for x > 0, the summands in Eq. (IA5l) decrease as i and j increase. 

We utilize the fact that, for a monotone decreasing continuous function F{x), a finite sum is 
bounded above by an integral as 

N „Ar 

^F{i)< F{x)dx + F{l). (A7) 
i=i 

Applying Eq. (IA7l) twice to Eq. (IA5l) and using pi = we have 

R'<^{ / G,i-—^x~^^y~^)dxdy + 2 (A8) 



The double integral in the bracket of Eq. dASI) is, by change of variables. 



I, = (X- l)2(iVe^-i)2 / / ^fA!^dudv, (A9) 



with A = 1 + and e = VKN'^/^'^ /(NifJ') ~ 0{N~^/^). Note that in Eq.(|A51) the upper limit 
of the integrals is eA^'' ~ 0(^]\[(3-^)/'^i>^-'^)^ and the front factor scales as 0{N^'^). We consider 
the three cases of A separately. 
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(i) When 2 < A < 3, since Gi{x) ~ x as x oo and ~ x^ as a; — > 0, the lower (upper) limit 
of the double integral in Eq. dAQI) can be expended to (oo) to give a finite value and hence 

/i < (A - ifiNe^-')' r r ^^^dudv ~ 0{N'-^). (AlO) 

Jo Jo {uv}^ 

(ii) When A = 3, the upper limit of the double integral is C(l) and the integrand near the lower 
limit behaves as ~ (nf )"^. We use < Gi{x) < for x > to get 

h < -{Ne^-^filnNf ~ 0{{\nN)^). (All) 
2 

(iii) When A > 3, proceeding as in the case of (ii), we find 

A<2(A-1)%^) -0{1). (A12) 
The single integral in the bracket of Eq . (IA8l) is, by change of variables, 

h = 2(A - 1)N5^~^ / -^ip^u, (A13) 



5 



with 5 = KN^~f'/('^{fi) ~ 0(iV^"^). Note that in Eq. dATSt the upper limit of the integrals is 
6N^^ ~ (9(j\r(3-^)/(->*-i)) j-j^g factor scales as 0{N'^~^). We proceed exactly the same as 
in the case of the double integral and find that 

(i) When 2 < A < 3, h < 2{\ - 1)N5^-^ ~ C>(Ar3-^). 

(ii) When A = 3, h < N6^-^ In ~ C(ln A^). 

(iii) When A > 3, h < j^^N6^ ~ C'(Ar"(A-3)/(A~i))_ 

Collecting these, we see that R' is bounded above as 

0{N^-^) if2<A<3, 
< ^ C((lnAr)2) if A = 3, (A14) 
C(l) ifA>3. 

Next we consider with n > 2. Similarly to Eq. (IA5l) . we have 
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Applying Eq. dXTt twice to Eq. (IA15l) . 

. f^.^ , NK 



Rn < 



r NK NK 

Vi Ji l^^^Q:^'^~'y~'^Tdxdy + 2j^ lGoi--^x-nrdx}, (A16) 



where Gq{x) = 1 — e~^. At this point, we use the piecewise linear upper bound for Gq{x) by 



Go 



X for < s < 1, 
1 for X > 1 . 



(A17) 



Since Gq{x) < Go{x) for x > 0, we can write Eq. (IA16l) as 



(uv) 



(A18) 



where e and 6 are defined above. Now the integrations in Eq. dAlSI) are elementary. Focusing only 
on the A^-dependences, we find that 

0{N^-^ In N) for2<A<3, 

0{{\nN)^) for A = Sand n = 2, 

C(l) for A = 3 and n > 3, 

C(A^2-n) for A > 3 and 2 < n < A - 1, 

C»((ln A^)2Ar2-") for A > 3 and n = A - 1, 

^ C(A^-"(^-3)/(A-i)) for A > 3 and n > A - 1. 

Putting these together, we finally have 



Rn. < < 



\R\<R' + y^< { 



n=2 



0{N^~^ In N) for2 < A < 3, 
0{{lnN)^) for A = 3, 
for A > 3. 



(A19) 



(A20) 



APPENDIX B: EVALUATION OF FINITE SUM IN GENERAL FORM 



In this APPENDIX B, we derive a general expansion formula for the sum 



1 ^ 

i=l 



(Bl) 



for small y{> 0) and ^ oo with Pi = i ^ / (n and A = l + l//i>2as before. We take F{x) 
to be a positive monotone increasing function which diverges slower than x^/'^ as x oo and has 
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an expansion F(x) = Yl'^=o fnx"- Converting the sum into an integral as in APPENDIX A, S{y) 
becomes, in the ^ oo limit, 

Siy) = (A - l)y^-' / ^dx. (B2) 

Jy X 

We first let A 7^ integer and ttiq < A < thq + 1 for some integer uiq. Then we define 

mo-l 

F{x) = F{x) - /n^" (B3) 

n=0 

and divide F{x) into two parts 

■mo-l 

F{x)= J2 /nx" + F(x), (B4) 

n=0 

Plugging Eq. (IB4l) into Eq. (IB2l) . the first finite sum can be integrated term by term to give 



S{y) = (A - 1) E Y^,y- + (A - 1).^-^ [ r ^dx - / 
Here we use the fact that F(x) ~ x*"" as a; ^ and hence 



x^ 



(B5) 



J(A) = (B6) 
Jo 

converges. The last term can now be integrated term by term using the expression of F. The result 



IS 



S{y) = (A - m^)y'-' - (A - 1) E n + l-\ y^^- ^^^^ 

n=0 

Note that F depends on mg, the integer part of A. When A = mo + 1 (integer), we set A = mo + 1— e 
in the above formula and let e — > 0+. In this way, the singular term obtains a logarithmic factor. 
The result is 

00 „ 

S{y) = m^iy^^ - m,f^,y^^ \ny - m, V —^2/", (B8) 

' ^ ri — nrrir. 

where 



n — mo 

n=0(7^mo) 



00 



, Fix) , Fix) - fmaX""" , 

^= / ^^dx+ / ^ ^ dx. (B9) 



A special case = 1 — exp(— x) has been treated 



inlbvl. 
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APPENDIX C: THE LEADING ORDER ANALYSIS OF Bn, 



,n2,n3,n4 



Bm,n2,ns,n4 is defined in Eq. (HI with M' = (A-2)M/(A-1) and = (A-2)Q/(A-1) and 
so on. To see how the leading order behavior of „2 ^,3 is determined, consider for simplicity 
the integral 



/oo 
dz z-^ tanh"! zM' tanW zQ' (CI) 

with the condition 1 > M' > Q' . 

When A is sufficiently large, the leading orders in M' and Q' are given by the first terms of the 
expansion of tanh x = x + • ■ ■ and we have 

/oo 
dz 

= a„,+„,M"^Q"^ (C2) 

Eq. (IC2l) with ai given in Eq. (l4TT) holds as long as A > rii + 722 + 1, but the integral in Eq. (IC2l) 
diverges when A < ni + 722 + 1 indicating appearance of the non-analytic term as the leading term. 

When A = rii + 722 + 1, the next leading order in Eq. (IC2l) cancels the divergence in a„i+„2 ^ 
give 

i3n„n,,o,o~ (A-l)(^)^^"'V"iQ"Mn(l/M). (C3) 
When 722 + 1 < A < ni + 722 + 1, one scales z — » zjM' in Eq. dCll) to find 



"OO 

-Bni.naAO = (A-l)M'^-i| / z"^ tanh"^ z tanh"^ (zQ/M) 

- y dz z-^ tanh'^i z tanh"^ {zQjM) | . (C4) 
The second term is (^(M^i+'^^+i-a^ smaller than the first whose leading contribution is 

/•oo 







C„„„2^^"'""'Q"', (C5) 



Cn,p= , . A-2 / c?a;x-^+Ptanh"a;, (C6) 







where C„ p, defined as 

(A -2)^-1 /-oo 
(A^ 

converges forp+l<A<r7, + p+ l. 
When A = 722 + 1, similarly to Eq. (IC3l) . 

'A-2\(^-i) 



Sn,n„o,o^(A-l)(-^) ' Q"Mn(M/Q). (C7) 
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When 1 < A < 722 + 1, one scales z ^ z/Q' m Eq. dCll) to write it as 

/>oo 

Bn,,n2flfl = (A-l)g'^"^ / d2z~^tanh"2 ztanh"^ (^Q/M). (C8) 

Jq 

Since Q <^ M <^ 1, ta,nh{zM/Q) ^ 1 for all z except near the origin where the contribution to 
the integral is negligible. Thus we have 

POO 

Bn„n,flfl ^ (A - l)Q'^-^ / dz z~^ tanh"^ z = Cn„oQ^~\ (C9) 

Jo 

The leading order terms of Bni,n2,m,n4 for various A's are listed in TABLE III. For simplicity, 
we do not show the A = integer cases in TABLE IIL For the border line cases of A dividing the 
regions of A with different expressions, a logarithm correction appears as given in Eq. (IB8I ) or (IC3I ) 
or (IC7I) . while for other integer values of A, the expressions are continuous. 
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integrals 


2 < A < 3 


3 < A < 4 


4 < A < 5 


5 < A < 6 


A > 6 


^2,1,0,0 


C2,iM^-2Q 


C2,iM^-2Q 


asM^Q 






^3,0,1,0 


C3,iM^-2Q3 




C3,iM^-2g3 






^2,2,0,0 




C2,2M^^3g2 


C2,2M^~3g2 






^0,3,0,0 












^4,0,0,1 






C4,iM^-2Q4 


C4,iM^-2Q4 




^0,2,0,1 












^0,4,0,0 













TABLE IH: The leading order terms of Eq.^. Here Cn,p = ((A - 2)^"V('^ - 
1)'*'^^) Jg°° dx x"'*'"'"^ tanh" x, which converges in the region ofp+1 < A < n + p + 1 and a; = 
(A-2)7[(A-l)'-i(A-l-0]. 
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